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Abstract 

We study the QCD color interactions between static two heavy quarks at zero temperature in 
a quenched SU (3) lattice gauge simulation: in addition to the standard singlet qq potentials, we 
calculate octet qq potentials, symmetric and antisymmetric qq potentials. It is shown that the 
antisymmetric qq channel behaves as a linearly rising potential at large quark separations. We 
further find that the qq octet and qq symmetric channels have the complex dependence on the 
distance; at short distances they are repulsive forces, while at large distances, they show linearly 
rising feature. Ratio of string tensions between qq singlet and qq antisymmetric potentials is 
described in terms of the Casimir factor. 
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I. INTRODUCTION 


Nonperturbative study of color QCD forces between static quarks is important for un¬ 
derstanding quark confinement and hadron phenomenology. The behavior of a heavy-quark 
potential that is defined by the closed gauge-invariant Wilson loop, has been studied ex¬ 


tensively in lattice simulations Q; the heavy-quark potential is a linearly rising potential, 
which can be explained with the Coulomb term and the linear term with the string tension. 
In most previous studies of heavy-quark potentials, a color singlet channel, which yields a 
physical potential, has been investigated. There are, however, several other color channels 
between two quarks. According to the £77(3) color decomposition for a quark-antiquark 
sector qq, 3 ® 3 = 1 © 8. The octet channel is significant for the extensive analysis of J/ij) 


photoproduction 


Furthermore, we have 3®3 = 3©6fora quark-quark sector qq , 


and in particular, the antisymmetric qq channel plays an essential role in the phenomenol¬ 


ogy of penta-quark hadrons 


is assumed jf|. Since the usual Wilson loop cannot give the color-decomposed potentials 
separately, the lattice study along this line was considered as very difficult. 


, in which the existence of a highly correlated diquark 


y cl: 

i. 


In order to understand quark-gluon plasma (QGP) physics |7j, the finite-temperature 
behavior of the QCD color forces between two static heavy quarks is studied in lattice 
simulations. At finite temperature, one expects that the linearly confining color force is 
weakened due to the screening effect in QGP. Quarks confined in hadrons are expected 
to move freely in the deconfinement phase after the QGP phase transition. This can be 
understood in the nonperturbative behavior of the screened heavy-quark potentials defined 
by a Polyakov line correlator, the form of which is reduced to a Yukawa-type potential with 
screening masses rather than the linearly rising potential. We therefore may expect the mass 


shi 


t and/or suppression of heavy quarkoniums as an indicator of the QGP phase transition 
9). However, recent heavy-ion collision experiments and the lattice QCD simulations 


13 


have indicated that this picture is too simple and QGP is a more complex system. 


Thus, in recent finite-temperature lattice simulations, the singlet and octet aa potentia 


WQ£\. 


and the symmetric and antisymmetric qq potentials have been investigated 
This is a progress comparing with previous lattice QCD studies, where the screened heavy- 
quark potential had focused on only a color-averaged channel. On the other hand, the 
nonperturbative behavior of the SU (3) color-decomposed channels at zero temperature has 


2 









not been investigated. 

Recently, using 577(2) lattice gauge simulation, Greensite, Olejnik and Zwanziger studied 
the color confinement mechanism at zero temperature and found that the Coulomb singlet 
heavy-quark potential in Coulomb gauge represents a strong linearly rising potential jlfil. (17], 
and this was later verified in SU(3) lattice gauge simulation Jl8(. The Coulomb singlet heavy- 
quark potential is calculated from partial-length Polyakov lines (PPL) in Coulomb gauge, 
not from the closed gauge-invariant Wilson loop. PPL is a Polyakov line with a restricted 


temporal extension 


lfl 
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19]. Thus, the combinations of the PPL correlators may enable 


us to carry out lattice calculations in the color-decomposed channels at zero temperature. 

In this letter, we study the behavior of the color-decomposed qq and qq potentials at 
zero temperature in the quenched SU(3) lattice gauge simulation. Applying the PPL cor¬ 
relator method Q, [l7| to the color-decomposed potentials, we carry out the 57/(3) lattice 
calculation of the octet, symmetric and antisymmetric potentials. We will present the first 
nonperturbative study of the antisymmetric qq potential at zero temperature. 

II. PARTIAL-LENGTH POLYAKOV LINE 


In this section, we give color decomposed correlators between two static heavy quarks 
and summarize how to fix the gauge on the lattice. 


A partial-length Polyakov line (PPL) can be defined as QQ 


L(x,T) =Y\_U 0 (x, t), T — 1,2, ■ ■ ■ , N t . 


( 1 ) 


4=1 


Here U 0 (x,t ) = exp(iagA 0 (x, t)) is a SU( 3) link variable in the temporal direction and a, 
g , A 0 (x,t ) and N t represent the lattice cutoff, the gauge coupling, the time component of 
gauge potential and the temporal-lattice size. A PPL correlator in the color-singlet channel 
is given by 

T) = j < Tr[L(R, T)L\ 0. T)]> , (2) 

where R stands for |x|. From Eq. 0 one evaluate a color-singlet potential on the lattice, 

Gi(R, T) 


V ^= l °Zl Gl (R,T + a)_ 

For the smallest temporal lattice extension, i.e., T = 0, we define 


( 3 ) 


V(R,0) = -\og[G 1 (R,l)]. 


( 4 ) 
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Greensite et al. argued that V(R, 0) in Coulomb gauge corresponds to a color-Coulomb 
potential lf| [l7|. The V(R,T) at T —> oo corresponds to a physical potential, which is 
usually calculated from the Wilson loops at T —> cx). These two potentials are expected to 
satisfy Zwanziger’s inequality, V p h ys (R ) < V cou i(R ) [21]. 

Now we extent this PPL correlator method to the other SU (3) color-decomposed poten¬ 


tials 


128 ]. A color-octet correlator for qq is given by 


G s (fl,T) = i<TrL(B,r)TrLt(0,T)>-l(rrL(fl,T)Lt(0,T)>, 


( 5 ) 


and qq correlators in the symmetric and antisymmetric channels are described as 


G e (R, T) = \ (TrL(R, T)TrL{ 0, T)) + ^ (TrL(R, T)L( 0, T)>, (6) 

Ga{R, T) = \ (TrL(R, T)TrL( 0, T)) - ^ (TrL(R, T)L{ 0, T )). (7) 

In the same way as described in Eqs. © and we will obtain the color-decomposed 
potentials in each color channel using the above correlators. The four potentials between two 
quarks are classified in terms of the Casimir factor on color SU (3) group in the fundamental 
representation: C qq = —| for a singlet, C qq = | for an octet, C qq = | for symmetric, 
Cq q = — | for antisymmetric; i.e., the qq singlet and qq antisymmetric channels yield an 
attractive force, whereas the qq octet and qq symmetric channels give a repulsive force. 

Since the color-decomposed potentials defined by the PPL correlators do not have gauge- 
invariant forms, we must fix the gauge. We use the Coulomb gauge realized on the lattice 
as 

3 

Max EE ReTr U}(x,t), (8) 

x i =1 

by repeating the following gauge rotations: 


Ui(x, t ) —> U“(x, t ) = u\x, t)Ui(x, t)u(x + i, t), 


( 9 ) 


where u G £77(3) is a gauge rotation matrix and U t (x, t) are spatial lattice link variables. 
Thus each lattice conhguration can be gauge fixed iteratively j22|. 

The Coulomb gauge fixing does not fix a gauge completely, and one can still perform a 
time-dependent gauge rotation on the Coulomb-gauge fixed links, 

Ui(x,t ) —> uj^(t)Ui(x,t)uj(t), 

Uo(x,t ) —> u^(t)Uo(x, t)uj[t + 1). 


( 10 ) 
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Thus, TrLTrZ/t and TrLL constructed by PPL, are not invariant under the transformation 
of Eq. (1101) . Therefore, when performing numerical simulations for the octet and two qq 
correlators including TrLTrJJ and Tr LL, we additionally implement a global temporal-gauge 
fixing on the Coulomb-gauge fixed configuration as 

Max ^ ^ ReTrf/p( x , t ) under Eq. (FTTH) . (11) 

X,t 

where V = N x N y N z is a spatial lattice volume and N x _ y>z is the spatial lattice size in each 
direction. Note that this gauge fixing does not affect intrinsic Coulomb gauge features. 


III. SIMULATION RESULTS 


We carry out 577(3) lattice gauge simulations in quench approximation to calculate the 
color decomposed PPL correlators. The lattice configurations are generated by the heat- 
bath Monte Carlo technique with a plaquette Wilson gauge action, and to fix a gauge we 


221 . 


adopt the iterative method 

Numerical results of the color decomposed potentials with T = 2 are displayed in Fig. |TJ 
This simulation is done on the 24 3 x 32 lattice at (3 — 5.90 (a ~ 0.124/m) and the 200 con¬ 
figurations measured by every 100 steps are used here. The singlet qq and the antisymmetric 
qq potentials rise linearly at large distances; i.e., they are an attractive confining potential. 
The octet qq and the symmetric qq channels show repulsive forces at short distances, while 
at large distances they have the linearly rising attractive force rather than the large-distance 


repulsive force. This nonperturbative behavior in 


the numerical result in 3d SU (2) gauge theory 


23 


he adjoint channel is comparable with 


24] 


In order to analyze the singlet qq and the antisymmetric qq potentials, we employ the 
following fitting function, 


V(R, T) = C + KR + A/R , A = —tt/12. 


( 12 ) 


where C is a constant and K corresponds to the string tension. This fitting works well: 
approximately y 2 /nd/ ~ 0(1) for the fitting range R= 2-6. The results are represented by 

solid lines in Fig. |T| 

We obtain the string tension Jig in the antisymmetric channel. The singlet-string tension 


K\ approaches asymptotically the Wilson-loop string tension 


18], but the T dependence 
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of Jig is unknown. Here, we investigate the ratio of the string tensions obtained from the 
singlet qq and the antisymmetric qq potentials; if the string tension itself is universal, then the 
difference between K\ and Jig can be understood in terms of the color SU(3) Casimir factor; 
that is, A'l/ilg —> C'A/C'| |J = (—4/3)/(—2/3) = 2. Those ratios calculated at (3 = 5.85 — 6.00 
are summarized in Fig. [21 The dependence of the ratio on /3 and T is weak and they are 
consistent with the expected values. We conclude that the relation between the singlet qq 
and the antisymmetric qq potentials is described in terms of the color 517(3) Casimir in 
nonperturbative regions. 

Since the repulsive channels have the complex dependence in the distance, we do not 
carry out a fitting analysis for them. 



FIG. 1: Singlet and octet potentials for qq sector and symmetric and antisymmetric potentials for 
qq sector. The solid curves represent the fitted results for the two attractive channels. The R with 
physical dimension is set as R = aR where the lattice cutoff a is approximately 0.124/m for this 
lattice simulation. 


IV. SUMMARY 

We have nonperturbatively studied the zero-temperature behavior of the color decom¬ 
posed potentials between static two heavy quarks in the quenched 517(3) lattice gauge 
simulation with the partial-length Polyakov line correlators. 
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Ratios, 24 x32, (3 =6.0 


Q--Q (i=6.00 
im| 3=5.90 
* «|3=5.95 
Ai|3=5.85 


FIG. 2: Ratios of string tensions in the singlet and antisymmetric channels. 


We show that the qq potential in the antisymmetric channel rises linearly at large dis¬ 
tances and is fitted by the function that includes the Coulomb and the string-tension linear 
terms. Moreover, the string tensions obtained from qq singlet and qq antisymmetric po¬ 
tentials seem to be understood in terms of the Casimir; in other words, the string tension 
itself is universal. As a result, not only the qq singlet channel but also the qq antisymmetric 
channel is found to be strongly correlated in nonperturbative regions. 

The octet qq and symmetric qq potentials show the complex behavior as a function of the 
distance in contrast to the case of the attractive channels: the qq octet and qq symmetric 
channels yield a repulsive force at short distances, while at large distances it seems to be 
reduced to a linear potential. The larger lattice simulation therefore may be required to deal 
with the long-range features in the repulsive channels. 

In the case of color non-singlet channels, it is especially important to investigate the 
volume size effect, because the color non-singlet potentials may include divergent parts in 
the infinite volume limit. 

Our numerical simulation has been concentrated on the two quarks system. With the 
PPL correlator we may deal with multiquark potentials on the lattice. The discovery of the 


penta-quark 


ladron 


state 


25, 


m 


I 


promotes the nonperturbative lattice research about the multiquark 


27]. The SU(3) lattice simulation of tetra- and penta-quark potentials 
also has been carried out H H 27] . The PPL correlator enables us to construct them 


without any assumption to make tentative forms of the multiquark potential, although the 
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gauge-fixing procedure is indispensable. 

The color decomposed potentials at zero temperature are computed in the quenched 
lattice simulation, and especially the behavior of the antisymmetric qq potential should 
be also investigated in a dynamical quark simulation. The PPL correlator ( or the usual 
Polyakov line correlator ) is applicable to the dynamical-quark lattice simulation. 
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